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Abstract 



We propose a simple model unifying two major approaches to the analysis 
of large multicomponent systems: interacting particle systems (IPS) and couple 
map lattices (CML) and show that in the weak interaction limit depending on 
fine properties of the interaction potential this model may demonstrate both con- 
densation/synchronization and independent motions. Note that one of the main 
paradigms of the CML theory is that the latter behavior is supposed to be generic. 
The model under consideration is related to dynamical networks and sheds a new 
light to the problem of synchronization under weak interactions. 

1 Introduction 

At present there are two major approaches to the analysis of large multicomponent sys- 
tems: Interacting Particle Systems (IPS) and Coupled Map Lattices (CML). Ideas and 
especially methods used in these approaches are strikingly different and it is hard to find 
a single article which discuss their connections on a reasonably serious level. Nevertheless 
in this paper we present a simple mathematical model enjoying interpretations in terms 
of both of these approaches which helps a lot in its study. 

Let A be a compact subset of M. d equipped with the Euclidian metric p(-, •) and let 
(T, A) := (g)f (T, A) be a direct product of N copies (in general < N < oo) of the mea- 
surable nonsingular dynamical system (T, A, B, fix)- Consider also a family of measurable 
maps Q £ : A — > A to which we shall refer as interactions. The parameter e > here 
measures how far the map Q e is from the identical map, namely e := sup sg ^p(Q £ !r, x), 
where p{x,y) := maxj p(x i: yi). 

Now a CML is defined as a composition of the direct product map and the interaction: 
T e := T o Q £ and the system without interactions can be formally written as T = T. 

A typical example of interactions considered in the literature is the so called "diffusive 
coupling" : 
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where the summation is taken over "neighboring elements" Jj of the local unit i. The 
formula (jl.ip indicates that speaking about a CML one usually has in mind a certain 
structure or topology of connections, say a fixed ordered graph of interactions (see e.g. 
[D El El [71 El [91 [ini [EU IE])- Another possibility is to consider a dynamic setup allowing 
a dynamical switching of connections. Some results in this direction were ob6ained in 
[T6] (in the case of time- varying couplings) and in |12j . In what follows we propose 
a simple model of a dynamical switching and show that in the weak interaction limit 
depending on fine properties of the interaction potential this model may demonstrate 
both condensation/synchronization and independent motions. 

Motivated by the idea of the interaction between the local systems induced by (some- 
what artificialjj "collisions" introduced in (12] we interpret the CML as follows. A con- 
figuration T e x := x* = {x\} at time t > is associated to the collection of identical 
particles with unit masses located at points x\ £ X which wander in the common domain 
X independently until they come close enough to each other. As an example one can 
think about a billiard type system. The interaction occurs only when the particles are 
coming e-close to each other and consists in the attraction to the local common center of 
gravity. 

Assume that (X, B,m, p) is a compact convex Lebesgue metric space X £ Mr. For a 
configuration x £ X and each index i denote by Jj = Ji(x) := {j : p(%i, Xj) < e} the set 
of indices of particles "interacting" with the i-th onej^l Then we define the z-th coordinate 
of the interaction map as follows: 



The parameter 7 £ [0, 1] and if 7 = only the second summand, corresponding the 
"center of gravity" of particles belonging to the ^-neighborhood of the point Xj, survives 
in the expression above. Therefore the physical meaning of the interaction under study 
is that the i-th particle is moved in the direction of the center of gravity of particles 
belonging to the e-neighborhood of the i-th particle. 

If 7 £ (0, 1) we say that the interaction is soft, 7 = corresponds to the situation 
when a particle jumps directly to the center of gravity, which we call rigid by contrast. 
Finally 7 = 1 nullifies the interaction. 

Our main results may be formulated as follows. Let the (T, X, B, /it) be a weakly 
mixing measurable dynamical system and let fiT be its only Sinai-Bowen-Ruelle (SBR) 
measure (see definitions and discussion in Section [2]). These assumptions are enough for 
the case of rigid interactions, but to study soft interactions one inevitably needs some 
additional smoothness type assumptions for the local map (see Section HJ). 

Theorem 1.1 For each N £ Z + there are constants et > 0, 7^ > such that V0<e< 
Et-i < 7 < 7t the CML (T e ,X) has a SBR measure Jt (which does not depend on e), 
is supported by the "diagonal" of X, and whose projections to each "coordinate" coincide 
with [At- If N = 2 this SBR measure is unique. 

1 The interaction discussed in [12] resembles more closely the exchange of velocities between colliding 
particles rather than the collision itself. The latter is modelled by a prior arrangement of non-overlapping 
"traps" for pairs of particles. 

2 Later we shall consider some other choices of the set Jj. 
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This result demonstrates that arbitrary weak interactions of a CML with a "generic" 
local map may demonstrate discontinuity at of the SBR measure /Z e considered as a 
function on e. Previously such results were known only for "wild" examples of local maps 
having periodic turning points type singularities, see e.g. [21 [j] . From the point of view of 
interacting particle systems this theorem can be interpreted as a kind of a condensation 
phenomenon when particles are gathering together under dynamics. 

In the opposite case when 7 is close to 1 we observe a more "classical" effect, namely 
under some reasonably general technical assumptions specified in Section [5] the only SBR 
measure of the CML converges weakly to the direct product measure as e — > 0. Bassicaly 
we need that the action of the transfer operator corresponding to the local map in a 
suitable Banach space of signed measures be quasi-compact, i.e. can be decomposed into 
a sum of a ^-contracting operator (with 9 < 1) and a finite dimensional projector. 

Theorem 1.2 For each N G Z + there is e' T > such that V < e < e' T and each 
7 G (2 N 6, 1] the CML (T e , X) has a unique SBR measure JI ea ^— ^> Jl T :— jjlt ® • • • ® l^r- 

N 

From the interacting particles interpretation this result tells that in this case weak 
interactions lead to independent particles motions. Thus soft interactions demonstrate 
very different and more "classical" effects in comparison to the rigid interactions. Note 
that here one may consider not only particles of different masses but even the local maps 
need not be identical. 

To understand better the nature of the phenomena under study let us consider the 
behavior of the CML in more detail and under a bit different point of view which is more 
closely related to the IPS approach. 

We say that a synchronization or condensation with the basin Y takes place in the 
CML (T E ,X) if limsup maxp(x*, x*) = for each configuration x G Y. Similarly a 

t— >oo M 

desynchronization with the basin Y means that liminf maxp(x',i') > for each x G Y. 

We say also that a certain property with the domain A and a probabilistic measure v is 
v- global if u(A) = 1. 

In words, the synchronization means that the dynamics converges to the diagonal 
D := {x G X : Xi = Xj Vz, j}, while the desynchronization corresponds to the absence of 
such convergence. We refer the reader to [15] for the discussion of various physical aspects 
of the synchronization. 

Let us start with the rigid case. The following statement formulated in terms of the 
de/synchronization not only clarifies the situation but also gives some additional insight 
about the dynamics. 

Theorem 1.3 Under the assumptions of Theorem li.il Vs > the synchronization with 
the open basin Y of positive product measure ~p T takes place. Moreover for N = 2 this 
synchronization is ~p, T - global. On the other hand, for N > 3 there is an analytic local map 
T for which the desynchronization with an open positive product measure ~p T basin takes 
place. 

It is worth note that our numerical simulations show that for all ID and 2D mixing 
maps that we tried for all initial particle configurations the global synchronization was 
observed. This is especially striking because first one might expect the synchronization 



only for almost all initial particle configurations, and second even for the cases where we 
are able to prove the desynchronization we do not see it in the numerical experiment. The 
explanation is that arbitrary small round-off errors may change the behavior of a chaotic 
dynamical system drastically (see e.g. PQ) and they are responsible in this case as well. 

The proof of Theorem 11.31 is based on the following observations. If two particles are 
coming e-close to each other, then after the interaction their positions coincide. There- 
fore if N = 2 we only need to show that for almost all 2-particle configurations their 
trajectories will hit simultaneously the e-neighborhood of the diagonal D. To prove the 
latter statement one uses that outside of this e-neighborhood there are no interactions and 
hence it is enough to check the same statement for the direct product system (T, X, B, /Z T ), 
which in turn follows from the weak mixing of the local map. 

One is tempted to extend this construction directly for the case N > 2 along the 
following lines. First, if the trajectory of an initial configuration hits an e-neighborhood 
of the diagonal D (as in the case N = 2) then after the interaction the synchronization 
takes place: all particles will share the same position. Second, additionally to the pre- 
vious argument one considers multiple (triple, etc.) collisions between particles (which 
occur near secondary diagonals of X) expecting that each collision reduces the number of 
unmatched particles. 

Unfortunately, there are two fundamental obstacles to this naive approach. First, a 
single collision of n > 2 particles does not necessarily imply that they will share the 
same location after the collision. Indeed, assume that we have n > 2 particles uniformly 
distributed along a circle of radius r n = -/ sin(7r/n). Then in the e- neighborhood of 
each particle there are two neighbors located at distance e. Therefore after the interac- 
tion instead of coming to a single common center of gravity, the particles will be again 
uniformly distributed along a circle of only a slightly smaller (provided n > 1) radius 
r' n := r(l + 2 cos(27r/n))/3, e.g. r' 4 = r 4 /3 and r' 5 rs r 5 x 0.539. Observe that we do 
not assume any smoothness of the local map T, and thus after the next application of T 
neighboring points might become arbitrary far from each other. 

The second and even more important obstacle is that even a single collision changes 
the trajectory of a particle and one cannot apply (at least directly) arguments related to 
ergodic properties of the original local system (T,X,B, fi). Moreover after each decrease 
of the number of unmatched particles one needs to check that the system is still in a 
"general position" in order to use again the mixing property. 

The result about the desynchronization shows that the smoothness of the map T does 
not cure these pathologies even in the simplest setting. In a sense the doubling map turns 
out to be the "worst" local one- dimensional map for our problem. 

Now we turn to the analysis of soft interactions. The most striking difference to 
the rigid case is that results of Theorem 11.11 or 11.31 are no longer available without some 
assumptions on the smoothness of the local system. Nevertheless Theorem 14. II (Section |4"]) 
shows that if the local map is Lipschits continuous then we recover results of condensation 
type when < 7 < 1. 

To study the opposite case when < 1 - 7 C 1 one needs to apply the transfer 
operators technics for suitable Banach spaces of signed measures. These matters will be 
discussed in Section using a multidimensional version of the Lasota-Yorke inequality. 
The idea here is to consider the action of interactions as weak perturbations to the dy- 
namics of each of local systems. This perturbative approach is well known but in our case 



there is an additional problem related to the fact that in distinction to the already studied 
situations the transfer operator related to the interaction might make a large contribution 
to the strong norm dpending on the total number of particles in the configuration. 

One might argue that if e = 1 and 7 1 then we recover the usual weakly intearcting 
CML setup. The assumption that e < 1 makes the already weak interactions to occur 
very rare and thus should preserve the "almost" direct structure of the invariant measure. 
Unfortunaly this "soft" argument does not work because the local structure of interactions 
leads to the creation of small regions where the "variation" of a measure might become 
arbitrary large under the action of interactions (see Section [5]). 

Additionally to the closeness in space the interaction might depend on the closeness in 
the "lattice position" . To be precise one considers the graph of interactions and assumes 
that it is locally finite, i.e. the degree of each vertex is finite. The degrees need not to be 
uniformly bounded. This allows to consider infinite systems, discussed briefly in Section [6] 
(see also a mean field approach in this Section). 

A few words about notation. We use a convention that for a Borel set A and a 
(signed/complex) measure \i the restriction of the measure \x to the set A is denoted by 
[a\a and /J,((p) :— J (p dp,. The bar notation x is used to mark variables describing the 
CML. Note also that | ■ | is used for very different objects throughout the paper and we 
follow the convention that in the case of a subset of integers \J\ means its cardinality, in 
the case of an interval - its length, in the case of a function \ip\ := esssup x and in 

the case of a (signed/complex) measure stands for its weak norm, while \\fi\\ stands 
for the strong norm (see Section loTTj) . 

2 Direct product systems: basic ergodic construc- 
tions 

Here we give a short description of standard definitions an constructions from ergodic 
theory which are necessary for the proof of our results. 

Recall that a measure n is T-invariant if and only if fj,(ip o T) — fjL(cp) for any \i- 
integrable function if : X — > M 1 . 

A measurable function if : X — > M 1 is called invariant with respect to a dynamical 
system (T, X, B, /1) (or simply T-invariant), if if = if o T almost everywhere with respect 
to the measure \i. 

A dynamical system (T, X, B, /1) is ergodic if each T-invariant function is a constant 
/x-a.e. 

A dynamical system (T, X, B, fi) is weak mixing if 

-. n—l 

- HT- k A nB)- n{A)n{B) I ™ V7L, B G B. 

A direct product of a pair of dynamical systems (T', X', £>', fi') and (T", X", B", /i") is 
a new dynamical system (T" ® T", X' £g> X", B' <g) B', // (g) //'), where the map T' ® T" : 
X' ® X" ->• X' ® X" is defined by the relation T ® T"(x',x") := (TV, TV), while all 
other objects are standard direct product of spaces, cr-algebras and measures respectively. 

By A N we denote the direct product of N identical sets A e B, and by (T® N , X N , B N , fx N ) 
- the direct product of N identical copies of a dynamical system (T,X,B, fi). 



Theorem 2.1 Let a dynamical system (T,X,B,fi) satisfy the weak mixing property. 
Then for any positive integer N < oo, measurable set A G B N with fi N (A) > and 
for almost any (with respect to the measure fi N ) collection x := {xi,X2, ■ ■ ■ ,%n} G X n 
there exists a moment of time t > such that T x G A. 

Proof. Let a dynamical system (r, Y, By, v) be ergodic. Then for any pair of measurable 
sets A,Be By with v{A)v{B) > there exists a positive integer k = k(A, B) < oo such 
that t r A n B 7^ 0. Indeed, assume that this is not true, i.e. r n A PI B — for any positive 
integer n. Consider a measurable set A^ := IJnez+ r?1 ^- Then ^(Aqc) > v(A) > 
and Aoo PI S = 0. Therefore the indicator function of a measurable set of positive 
i/-measure is r-invariant but is not a constant a.e. which contradicts to the ergodicity. 

Therefore it is enough to show that the dynamical system (T® , X , B , fi ) is er- 
godic. For that we shall take advantage of the fact that the weak mixing property is 
preserved under the direct product of weak mixing dynamical systems (see e.g. [12]). To 
complete the proof its remains to note that the weak mixing implies ergodicity. □ 

It is of interest that one cannot weaken the conditions of Theorem 12.11 replacing the 
weak mixing of the original dynamical system by the ergodicity. The problem is that the 
direct product of ergodic dynamical systems needs not to be ergodic as well: consider a 
direct product of two identical irrational unit circle rotations. On the other hand, the 
weak mixing condition is not necessary as well (see |3] for an example of a nonergodic 
dynamical system for which every open neighborhood A of the diagonal of X x X satisfies 
the claim of Theorem 12. ip . 

In what follows we often deal with the e-neighborhood of the diagonal D of X = X N 
which we denote by D £ . It is not difficult to calculate the Lebesgue measure of this set 
but the product measure ~p(D e ) may vary sensitively with /i. To estimate it we use the 
techniques of measurable partitions. 

Recall that a measurable partition of (X,B) is a collection A := {Aj}, Aj e B such 
that Aj n Aj = 0, flj Aj = X. The diameter of a partition is the largest diameter of its 
elements. 

Lemma 2.1 For each < e and a probabilistic measure /i we have < ~p{D £ ) < 1, 
moreover sup M /Z(D £ ) = 1. Assume now that > there exists an partition of X of 
diameter e with cardinality n e < Ce~ d . Then inf^/Z {D E ) > {Ce~ d )~ N+1 . 

Proof. Let A be a finite partition of X of diameter e > 0. Since X is compact such 
partitions always exist. Denote aj := /i(Aj) then J2i a i = M^O = 1- Observe now that 
Af^ := Aj x . . . x Aj C D £ . Therefore the lower estimate oifl{D e ) follows from the trivial 

N_ 

inequality Jl(D £ ) > EM^i)) N = E^f > 0. 

Now let /i be concentrated at a single point x G X, i.e. /i({x}) = 1. Then obviously 
■p(D e ) e 1 Ve > 0. 

It remains to prove the lower estimate under the additional assumption about the 

cardinality of the partition A. Consider a power average SV({ai}) := \^-YH=i a f ) 
of n £ nonnegative entries aj. It is known that Sn(-) is monotonous on the nonnegative 



parameter N. Thus Sn({ch}) > S\({ai}) = — . Therefore n(D e ) > 'Yl li {^{^-i)) N = 
£,af = n e (S N ({ ai })) N > n e (S 1 {{a i })) N = n~A ~ ' □ 



Lemma 2.2 Let A C X with Jl(A) = 0. Then Jl((T)~ l A) + Jl((Q £ )~ l A) = Ve,7 > 0. 

Proof. The fact that Ji({T)~ l A) = follows from the nonsingularity of the local map 
T. If e — the interaction does not occur and hence Q is an identical map. Therefore 
it is enough to consider e > 0. ViV G Z + and Va; G X the interaction map Q e may be 
written as a nonsingular linear map. However the matrix defining this map depends on 
the "grouping" Ji(x). Nevertheless for the finite number N of particles the total number 
of various "grouping" is finite and hence Q £ may have only a finite number of nonsingular 
representations. □ 

A map T : A — » A induces the transfer operator T* acting in the space of signed 
measures (generalized functions) Ai on X by the formula T*fi(A) := fi(T~ 1 A) for each 
A G B and /i G M. From this point of view a measure \x is T- invariant if and only if 
T*H = (i. 

A probabilistic measure G M. is called the Sinai-Bowen-Ruelle (SBR) measure 
for the dynamical system (T, X, B) and a reference measure m (say Lebesgue measure 
on X) if there exists an open subset Y C A such that for each probabilistic measure 
/i G absolutely continuous with respect to m and such that f-i(Y) = 1 we have weak 
convergence - Y^=o T* 1 fi n —^¥ \xt- The set Y is called the frasm of attraction for the 
measure //y Obviously an SBR measure is T- invariant. 

There are also different approaches for the definition of the SBR measure and we refer 
the reader to [5] for their discussion and conditions under which those approaches agree 
with each other. 



3 Rigid interactions 
3.1 Proof of Theorem IT51 

Let N = 2. Consider the e- neighborhood of the diagonal D £ in A. By Lemma 12.11 
~p(D e ) > 0. Therefore by Theorem 12.11 for G A there exists the first moment of 

_ t(x)_ 

time < t(x) < oo such that T x G D £ . Denote the set of full /i-measure for which this 
holds by X £ . 

By the definition of the rigid interaction, Q £ D £ C D while Q £ x = x if x ^ D e . 

Therefore Vx G A e , t G {0, . . . ,t(x) — 1} we have T £ x = T*x and T^x G -D which 
proves the global synchronization if N = 2. 

Local synchronization for an arbitrary N > 2 follows from the invariance of the e- 
neighborhood of the diagonal with respect to the dynamics. However, if N > 3 the 
observation that T £ x = T x V£ G {0, . . . , t(x) — 1} does not hold. 

It remains to show that when N > 3 even the analyticity of the map T does not 
guarantee the global synchronization. 



Lemma 3.1 Let X := S 1 (unit circle), the local system be governed by the doubling map 
Tx := {2x} and let N = 3. Then V < £ <C 1 the desynchronization with the domain of 
positive Lebesgue measure occurs. 

Proof. Let ij 6 I, i G {1,2,3} and denote a := x 2 — £i, b := x 3 — x 2 . Denote by 
A the subset of X for which < a, b < e and a + b > e. Under the interaction the 
coordinates X{ will be changed to x\ such that the distances between them will become 
equal to a! := (a + 2b)/6 and b' : = (2a + b)/6. Applying the doubling map we get the 
new pair of distances a" := (a + 26)/3 and b" := (2a + 6)/3. Since < a", 6" < e and 
a" + b" = a + b > e the new configuration again belongs to the set A. Now the observation 
that the product Lebesgue measure m(A) = e 2 /2 > finishes the proof. □ 

The extension of Lemma [3. II for the case when dim(X) > 1 is straightforward. More- 
over, the discussion after the formulation of Theorem 11.31 demonstrates that for large 
TV the local map demonstrating the desynchronization may be chosen o(l/iV)-close to 
identical for each e > 0. 

3.2 Proof of Theorem LTTT1 

One might think that Theorem 11.11 is a direct consequence of Theorem 11.31 Indeed by 
Theorem 11.31 a.a. trajectories of our CML Vfe > after a finite number of iterations hit 
the diagonal. Hence any forward invariant set of T e belongs to the diagonal D and obvi- 
ously the analysis of the invariant measures may be restricted to the "one-dimensional" 
dynamics on the diagonal. On the other hand, we do not assume that the T-invariant 
measure /i is unique and hence typically there is a subset Y C X on untypical points 
leading to statistics different from /i. Still the //-measure of this exceptional set is zero. 
Therefore using that number of iterations before to hit the diagonal is finite (but not 
uniformly bounded) and the result of Lemma 12.21 we deduce that the set of ~p- "typical" 
points is of full /I-measure. 

3.3 The closest rigid interactions 

As usual only the particles from the e- neighborhood of Xi will be included to Ji(x), but 
now we consider a special (simplified) choice of 

Ji(x) := {j : p(xi,Xj) < e, p{x h Xj) = min p(x h x k )}, 

i.e. this collection contains only the closest particles to the i-th one. We shall refer to this 
choice of Jiix) as the closest interaction. Note that the set of configurations x for which 
maxj |</j(a?)| > 2 has Lebesgue measure 0. 

Theorem 3.1 Let the (T,X,B,fi) be a weakly mixing measurable DS and /x be its only 
SBR measure. Assume also that the interaction is rigid with the above choice of Ji(x). 
Then for each N e Z + , e > the JZ-global synchronization takes place. 

Proof. For a configuration x G X denote by £(x) the minimum distance between 
particles in x. Observe that for /Z-a.a. configurations each Ji(x) consists of at most one 



particle and the minimal distance is achieved at a single pair of particles. By the definition 
of the interaction if £(x) < e these two particles may interact only between themselves 
and hence their positions after the interaction will coincide with the common center of 
gravity. 

On the other hand, if £(x) > e no interactions occur and one may use the same 
argument as in the proof of Theorem 11.31 to show that /Z-a.s. this effect takes place. 

To finalize the inductive construction, observe that for each particle the matching 
with some other particles takes place after a finite number of time steps and hence by 
Lemma 12.21 the Jt measure of "non generic" initial configurations leading to the non- 
uniqueness of £(x) is zero. □ 



4 Soft interactions. Case < 7 <C 1 

In this Section we study the intermediate case when < 7 < 1. It is easy to see that if 
7 is close enough to 1 then we are basically in the same situation as in the case of rigid 
interactions. 

Theorem 4.1 Let the (T,X,B,p) be a weakly mixing measurable DS and p be its only 
SBR measure. Assume also that there is a constant < A < I/7 such that p(Tx,Ty) < 
Ap(x,y) Vx, y G X Then the claims of Theorems \l.l\ and \1.3\ hold true. 

Proof. Let x G D £ . Then p(xi,Xj) < e V2, j and hence all "particles" do interact with 
each other. Denote by z their common center of gravity. Then (Q £ x)i := jXi + (1 — j)z 
which implies that 

p((QsX)h (Qs~x)j) < lp(xi,Xj) 

Thus Q £ D £ CD 7 , 

On the other hand TD £ C D?y £ by the assumption on the map T. Therefore 

ToQ £ D £ CD Al£ CD £ . 

Moreover, 

p«? e x) h (Tlx),) < (A 7 )V(zi,*i) ^ Wi,j. 

The completion of the proof follows exactly to the same arguments as in the case of 
rigid interactions. □ 

Corollary 4.1 Let At be the modulus of the largest Lyapunov multiplier of the map T. 
Then the claims of Theorems \1.1\ and \1.3\ hold true if and only if < 70 := l/A-p. 

Proof. The direct statement follows from the argument above applied to D £ /\ T rather 
than to D £ . To prove the inverse statement one observes that the modulus of the largest 
Lyapunov multiplier of the map T o Q £ cannot be smaller than A7-7 > 1. □ 

Remark. Despite the similarity between the case under consideration and the case of 
rigid interactions there is an important difference in that once x G D £ all "particles" 
will immediately form a cluster on the next time step in the rigid case, while an infinite 
number of iterations is nessasary for this in the soft case. 



5 Soft interactions. Case < 1 — 7 1 



Since 7 = 1 correspond to the absence of interactions (and thus the SBR measure of the 
multicomponent system is equal to the direct product of local SBR measures) one expects 
to observe a kind of phase transition when the parameter 7 grows from to 1. In what 
follows we shall study what happens when the parameter 7 becomes very close to 1. 

In the Introduction we already mentioned that for our purpose it is enough to assume 
that the local transfer operator T* be quasi-compact in a suitable Banach space of signed 
measures. This means that T* may be represented as a sum of a contraction and a 
compact operator. Below we shall show that this property is satisfied e.g. for the so 
called piecewise-expanding maps. 

To this end we need to introduce a proper Banach space of signed measures and to 
describe their properties. 

5.1 Transfer-operator approach and BV measures 

Recall that X is a unit iV-dimensional Eucledian cube equipped with the standard Borel a- 
algebra B = B N . The map T : X — > X induces the transfer operator T acting in the space 
of signed measures (generalized functions) Ai on X by the formula T* jj(A) := jj,(T A) 
for each A e B and /Z 6 M. 

Starting from [HI E] the approach to the analysis of transfer operators in terms of the 
so called dual norms proved to be efficient and became popular. To introduce the dual 
norms in the space of signed measures (generalized functions) on X we start with spaces 
of test-functions 

T := {if E C\X) : \ip\ < 1}, Jo := {<? & C\X) : \ip\ < 1, <p\ _ = 0}, 
T L := {<£ G L(X) : ^ < 1}. 
The following are two versions of the "variation" of a signed measure: 
V(/l) := max sup ~p(diip) , V (/l) := max sup Ji(di(p). 

The latter functional gives the variation of the density of the measure \i with respect to the 
Lebesgue measure, while the former also gives the variation of this density but considered 
as a function from ~R N taking zero value outside of X. An important advantage of the 
above definition of the variation is that in the case of a measure having the direct product 
structure properties of its variation can be easily obtained from their one-dimensional 
counterparts. Therefore we shall give proofs only for one- dimensional statements and 
refer the reader e.g. to [H] for the multidimensional setting. 
Define also the L 1 -norm of the (signed) measure: 

|/Z| := sup Jl(<p) 

which we shall call the weak norm. 

Lemma 5.1 (a) V(/Z|y) < V(/Z); 

(b) Vo(/Z| F ) < V(7Z |F ) < 2V (/Z| F ) + 2\jI lY \/m{Y); 

(c) if Y is a proper rectangl^ then < | m(Y) V(/Z), in particular, \p\ < |V(/Z). 
3 i.e. a direct product of n intervals. 



Proof, (a) V(/i|y) = sup veT n\ Y (p) < sup veF fx(p) = V(/i). 

(b) For if> G T set <^ := x((p(x) — y(l)) + (1 — x)(<^(x) — y(0)). Then |<^ | < 2|<p| 
and since <^o(0) = (fo(l) = we obtain that hcpo G Tq. Therefore fJ>((p') = Mv^o + (y?(l) ~~ 

< 2Vo(/x) + 2|/i|, which proves the inequality for the case Y — X. The general 
case can be proven similarly. 

(c) Decompose the signed measure fi\y '■— /■*+ — fJ>- into positive and negative compo- 
nents and set Y± := supp(/^±). The function <p(x) := m(Y + n[0,x])—m(Y_n[0,x]) — ^m(Y) 
is continuous on X and \<p\ < ^m(Y). On the other hand, by definition \fj,\y\ = My 9 ') — 
|m(y)V(/x) since 2(p/m(Y) G J 7 is a valid test-function. □ 



Therefore the functional V(/l) is actually a norm (which we denote by ||/l||) and is 
equivalent to a more common strong norm V (/x) + \fi\- Therefore we shall refer to 1 1 • 1 1 as 
a strong norm. Note also that for the Lebesgue measure V(m) = 2. The set of (signed) 
measures /i with < oo we shall call measures of bounded variation and denote this 
set by BV. 

Using this terminology we may rewrite Theorem 11.21 claiming the convergence to the 
direct product measure as follows. 

Theorem 5.1 Let the map T have the only one SBR measure ht, and let there are 
constants 0<#<l<O<oo such that 

\\T*fi\\ < 0|H + @H (5-1) 

for each fi GBV. If2 N 9/^ < 1 then for each < e < 1 the CML (T E ,X) has the only 
one SBR measure ~Ji e ^> JI T - the direct product of the local SBR measures. 

To give a specific model satisfying to our assumptions consider the class of piecewise 
expanding maps. Let X : = [0, 1] and r : [0, 1] —¥ [0, 1] be a piecewise C 2 -smooth map, i.e. 
there is a finite partition of X into intervals X{ on each of which the map r is bijective, 
C 2 -smooth, and inf a . \r'(x)\ > A > 0, /?i(r) := Amin 2 . /3 2 (t) := sup x \(l/r'(x))'\ < oo. 
Such maps are called X- expanding. Set /3(r) := A(t) + /M 7 ")- 

Lemma 5.2 (Lasota-Yorke inequality) Let the maps r 1; t 2 , . . . , be X-expanding, r stands 
for their direct product, and let /3(r) := maxj /3(Tj). T/ien 

||r*7Z|| < |||7I|| +/3(T) \% (5.2) 

Proof. First observe that from Lemma [5.11 (b) it follows that for each ip G T 



li(<f/) < 2\(p\ (VM + |/i 



(5.3) 



As (<p o r)'(ic) = <£>'(t(;e)) ■ t'(x) for each \ (UjSXj), we have 

r*/i(¥0 = r) = ° r)'A0 = M«V? ° r)/r)') - n({<p o r) ■ (1/r)'). 

To estimate the first term we apply (15.31) . while the second term is bounded by fair) 
□ 




Figure 1: Local structure of the map r (left). Decomposition of the test-function (right). 



Corollary 5.3 Under the assumptions of Lemma \5.S\ there exists a probabilistic r -invariant 
measure 

Proof. Choose k G Z + large enough such that \ k > 2. Denote r := 2/X k < 1 and let 
~p G BV be a probabilistic measure. Then for each n G Z + we have: 

II — *nk / — \|| ^ ran — II , P\T) \ — | 

1 — r 

Thus the sequence ji n := T* nk JI satisfies the conditions of the embedding of BV to the 
set of measures having absolutely continuous densities with respect to Lebesgue measure, 
which we denote by L 1 , and, hence, there exists the limit point of this sequence, i.e. 
JI n . /Zqo G BV with H/iooll < yr^|/i|- O n the other hand, being a limit point of this 
sequence the measure /Zoo satisfies the relation r*/!^ = /Z^ an d thus is r-invariant. □ 



Corollary 5.4 The transfer operator corresponding to the map T under an additional 
assumption of the uniqueness of SBR measures for each map T{ satisfies the conditions of 
Theorem \5.1\ with 6 := 2/A and 9 := /3(r). 

Consider now a more specific family of piecewise C 2 -smooth maps r : X — > X. For a 
given positive integer n let < a\ < a[ < < a' 2 < . . . < a n < a' n < 1 and set A { := 
[di, a-], i G {1, 2, ... , n}, A := UjAj, and B := X \ A. These intervals define a partition 

{T' CC ' — ■ £C | (2, ■ if £ ■ 
1 % . % , where r,A C A, a» > Vi. 

x otherwise 

Fig. [1] (left) demonstrates the shape of r in a neighborhood of an interval Ai. 

We shall be interested in the properties of the transfer-operator r* when \A\ 1. Due 

to this restriction the application of Theorem 15.21 gives an estimate with the second term 

going to infinity as the diameter of the partition goes to zero. In order to overcome this 

difficulty we develop a new approach to estimate the norm of the transfer operator in this 

case. 

Lemma 5.5 ||r>|| < (n+ 1 + Ei — ) and \fi-r*fi\ < H n + 2 + J2i ~) m ( A ) INI- 



Proof. The idea used in the proof of Lemma [5711 (b) is to interpolate linearly between the 
values of the test-function at boundary points of the partition {X{\ and to estimate the 
contribution of this interpolation into the integral against the weak norm, rather than the 
strong one. In the case under consideration the lengths of the intervals of monotonicity 
might be arbitrary small which does not allow to apply his trick. Instead we shall treat 
each interval of monotonicity separately extending the test-function by two constants 
equal to the values at boundary points outside of the interval (see Fig. [1] (right)). 

Observe that for each i the function Ti(x) can be extended as a linear function to the 
whole X. For a test-function ip G J 7 and a (signed) measure /i we have 

T*n(<ft) = fi((p' o r) = fi\ B {(p') + H\Ai(<p' ° t) 

i 

= P\b(<p') + H\a(<p') ~ P\a(<p') + Yl WAiW r ) 

i 

n 

= riff) + Y {WAiW ° T ) ~ P\Ai(<P')) ■ 
i 

Let x G Ai then </?' o r = (</? or)' • (r,-) -1 = ^-((/jo t)' and ip o r G T is a valid test-function. 
Therefore ^aX 1 ^' or ) — Summing up all contributions from the integrations over 

Ai we get 



*W)<(i+»+Ei) 



l/'l 



since //^(vO < I HI Vz by Lemma |5~T1 

To estimate |/x — r*/x| observe that the measures differ only on the intervals Ai. Set 
u := fj, — r*fi and A := UjAj. Then 

M = M < ^m(A)V(i/) < ±m(A)(|M| + ||r>||) 



< -mi A) \ n + 2 + V— | 

" 2 v r<v 



□ 



It might seem that the multiplier (n + 1 + — ) is overpessimistic, but the trivial 
example of the Lebesgue measure m on X immediately shows that ||r*m|| = 2(n + 1 + 

E. - .) while llmll = 2. 

As we shall see the argument used in this proof is the key point in the proof of 
Theorem 15.11 Apart from this Lemma 15.51 allows to apply the operator approach to a new 
class of small but discontinuous perturbations. 

Theorem 5.2 Let T : X — > X be a X-expanding map with A > 1 having a unique SBR 
measure and let r be a piecewise linear map described above with \A\ = 5 and a given 
collection of slopes {a{\ such that 2{n + 1 + ^ — ) < A. Then for each < 5 <C 1 the 

dynamical system (r o T,X) has a unique smooth invariant measure fxs A*t- 



Proof. Combining results of Lemmas 15.21 15.51 we get 



^ OLi 

% 

s(»+i+E^)5iMi + (»+i+E^)/'WW- 

Therefore we are in a position to apply Corollary 15.31 to the map r o T which yields the 
existence of the probabilistic invariant measure fj,s G BV. 
On the other hand, by Lemma 15.51 

\ji — r*fi\ < C6\\(j,\ | 

which by the now standard perturbation argument (see e.g. [HE]) implies the convergence 
Us — > ht- □ 

To finish this preparatory part let us formulate an estimate of the action of the transfer 
operator of a contracting afline map which can be proven by a direct inspection. 

Lemma 5.6 Let X be the N -dimensional unit cube and letr(x) := Gx + H be an affine 
map from X into itself such that £{G^) > ot£(£) for each £ G M. N and any norm £(■). Then 
||r*7Z||<l/a||7l||. 

5.2 Proof of Theorem 15.11 

One is tempted to argue as follows. Assume that the parameter e is of order of the 
diameter of the set X. Then each pair of particles is interacting between themselves and 
we are coming to the well known mean field interaction model. One can show that when 
7 goes to 1 the only SBR measure of the mean field model converges to the direct product 
measure, i.e. the subsystems behave independently. Now the decrease of e leads only to 
the decrease of the frequency of interactions between particles. Thus already "almost" 
independent particle are becoming even more independent. In fact the situation is much 
more complicated. The point is that the existence of a large number of small islands in the 
phase space where different combinatorial types of interactions actually take place leads 
to a severe amplification of the variation of a measure under the action of the operator 
Q e . Consider this in detail. 

The definition of the map Q e implies that configurations x G X whose coordinates 
have pair distances larger than e are fixed points of the map Q e . On the remaining part 
of the phase space X consisting of a large (of order 2^) number of disjoint components 
of small Lebesgue measure the map Q £ is linear and contracting in each of them. Thus 
the structure of the multi-dimensional map Q e is very similar to the one-dimensional map 
t considered in Lemma 15.51 Therefore we shall use basically the same strategy for the 
proof. 

For a given test-function (f G T we need to evaluate the functional maxj Qji[di^p). 
Fix some index i and for each subset J of different integers belonging to the set 
{1,2,..., iV} and containing the index i define a set 

Aj := {x G X : \x,i — xA < e Vj G J, [x, — x k \ > e Vfc ^ J}, 



and let A := Utj^xAj, B := X \ A. Then the interaction with the i-th particle occurs 
only for I6i and the sets A, B define a finite partition of X. 
We have 

Ql~P(di<f) = //(%> O Q £ ) = Ti\ B (di<p °Q £ )+Y1 hAj&V ° Qe) 

\j\>i 

\J\>1 

= JJ{di(f) + Y (p\Aj{dM ° Qe) - Ji\Aj(di<p)) 
\J\>1 

Denote by <p' the vector of partial derivatives of <p, by Q £ the matrix of partial derivatives 
of the map Q £ , and by (g^) the matrix inverse to the matrix Q e (i.e. (g^) = (Q £ ) _1 ). 
Then 

(dicp) o Q e = {(cp o Q £ ) ■ (Q' £ )- l )i = ^ dj(<P ° Qe) ■ Qji 

j 

= Yl d 3^P ° Qe ■ Qji) ~ ? ° Qz ' d i^- 
3 3 

Denoting ip^ := ip o Q ■ e C 1 we rewrite the expression for the action of the transfer- 
operator on a measure restricted to Aj as follows: 

~P\Aj( d i<P Qe) < $^ sup ' V (/ J ) + JZ sup ' \"P\- 

X . X 

3 3 

The interaction inside of each region Aj is described by a linear function and thus the 
last term is equal to zero for x G Aj. 

In general the prefactor V. sup x |<^i(a;)| can be estimated as 

^sup|^i(x)| < ^sup|((g^)~ 1 ) ii |. 

X . X 

3 3 

In our case the upper estimate can be done explicitly. Observe that for a given set Aj 
the restriction of the map Q e to Aj is a linear map (which we denote by L) has a very 
simple structure. Namely, considering only "interacting coordinates" one can rewrite this 
map as an affine map Ly := Gy + H with G := 7 + l^-E. Here n stands for the number 
of the "interacting coordinates" and all entries of the n x n matrix E are ones. It is easy 
to check that £(G£) > j£(£) for any norm £(■) and (6l" and the equality is achieved 
on a vector £ having the only one nontrivial coordinate. Therefore using Lemma 15.61 and 
setting a := 7 we get ||£*/Z|| < 1/7- Thus the prefactor can be estimated from above as 
I/7 uniformly on J. 

Now using that the number of different collections J cannot exceed 2^ we get 

I, 2 9.. I, 

WQeVl < \W\- 

7 

Combining this result with the Lasota-Yorke type inequality (15.1 p for the direct product 
map we estimate the strong norm of the action of the transfer-operator of the CML as 

\\tji\\ < e\\Q*ji\\ + e\Q* £ ji\ < — ||/i|| + e\jt\. 



Now using again the same trick as in the proof of Lemma 15.51 we show that the measures 
Tjl and T Ji are close in the weak (L 1 ) norm for small enough e > 0. 
For a collection of indices / introduce a set 

Bj := {x G X : Vii G I 3i 2 , . . . , i k G I : \x ii - x ii+1 \ < e, min |xj - xA > e}. 

In words, the configurations from the set Bj satisfy the condition that all /-particles (i.e. 
those with indices from J) are connected by e-chains, while all others are far enough. 

The sets {Bj} define a finite partition of X and the map Q e differs from the identical 
map only on the sets Bj with |/| > 1. Thus the signed measure V := Qjl — Jl is supported 
only on the sets Bj with |/| > 1. Denote Bj := {x G X : \xi — Xj\ < e\I\ Vi,j G /}. 
Obviously Bj is a proper rectangle and Bj C Bj V/. Applying Lemma [5.1( c) we get 

M = I^IUifl^Brl = ^2W\B I \<Y1 \ V \B T \ 
\I\>1 |/|>1 

IJI>1 

<-\T. w\r (i + ^)i W i^o. 

IJI>1 v ' 7 

The completion of the proof follows the same perturbation argument as in the proof 
of Theorem 15.21 □ 



6 Generalizations 

1. General interaction potential. So far we have considered local interactions based 
on a somewhat non-physical model of attraction. Indeed, normally by the attraction one 
means something more close to the gravitation law. In order to include these more general 
local interactions (excluded by (11. 2p ) consider an "interaction potential" U : X — Y M and 
set U e (x) := U(x/e). Then one defines the following generalization of the dynamically 
switching interactions: 

1 — 7 

{Q £ x)i := jxi + -rrr2_^ x 3 ' U ^ x i " ( 6 - 1 ) 

The interaction in (11.21) corresponds to the potential U(x) defined by the indicator function 
1 [—1,1] of the interval [—1, 1]. Assuming that U(x) > and making some regularity type 
assumptions on the potential one recovers all results obtained in Theorems 11.11 and 11.21 

2. Random local dynamics. It worth notice that the dynamics of the local units of 
the multicomponent system under consideration needs not to be deterministic. Indeed, it 
might be defined by a stochastic Markov chain satisfying the weak mixing condition (in 
Theorem 11.11) and some additional assumptions about the induced operator acting in the 
space of signed measures (in Theorem 11.21) . 



3. Infinite particle systems. Strictly speaking our definition of the interaction does 
not allow to consider infinite particle systems since the notion of the center of gravity is 
not well defined in this case. To overcome this difficulty one may assume that there exists 
a certain locally finite graph of interactions and that the dynamical switching occurs 
only between neighboring elements in this graph. This means that the sets J; satisfy the 
property that only elements J connected to i may belong to them. All our results can be 
extended to this setup. 

Another and potentially more promising approach (at least in the "independent phase" ) 
is to consider a mean field approximation scheme. Let /i be a probability distribution de- 
scribing the position of a particle (and assume that it is the same for all particles). If a 
given particle is located at a point x G X then the mean field approximation allows to 
calculate the center of gravity of the particles in the e-neighborhood B e (x) of this point 
as (i(b~ (x)) Ib Myd^iv)- Therefore one rewrites the interaction as 

1-7 f 

Qe a ^ x ■= i x + ,t> , ^ / ydfi(y)- 



fi(B e (x)) J B 



Denoting (for a given /x) by Q* i7iM the induced action of the map Q £ , 7lAt in the space of 
signed measures we obtain the description of the mean field approximation in this space: 

T e > := Q:, 7iT> T>. (6.2) 

In distinction to the transfer operators considered in the previous Sections the operator 
T e * 7 is nonlinear which complicates its analysis a lot. 

In the simplest case when X := S 1 and Tx := 2x it is easy to show that the Lebesgue 
measure is T e * 7 -invariant for all e, 7 and that any measure uniformly distributed on a 
periodic trajectory is invariant for small enough e > 0. Nevertheless the analysis of 
stability of these measures (i.e. the construction of the analogue of the SBR measure) is 
a much more delicate task. Even in this simple example one needs to develop a special 
technique to study properties of the nonlinear transfer operator. Therefore this analysis 
will be discussed in a separate publication 
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